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1  INTRODUCTION 

Exploitation  of  the  phenomena  of  acousto-optic  diffraction  to  construct  various  real-time 
signal  processing  architectures  is  a  familiar  principle  and  a  well  established  technology. 

The  ultimate  microwave  performance  of  these  devices  (in  terms  of  achieving  gigahertz 
operating  frequencies,  large  bandwidths  and  high  resolution)  has  been  severely  limited  by 
exponential  attenuation  of  the  signal  due  to  scattering  of  the  coherent  signal  phonons  by 
incoherent  thermal  phonons  in  the  acousto-optic  material.  The  room  temperature  co^ 
frequency  dependence  of  the  attenuation  coefficient  typically  restricts  usable  performance  to  a 
maximum  frequency  of  less  than  3  GHz,  with  a  1.2  GHz  bandwidth. 

Recent  research  into  the  characteristics  of  cryogenically  cooled  Bragg  deflectors  [1,  2]  has 
shown  that  efficient  operation  of  these  devices  at  frequencies  well  into  the  gigahertz  regime 
is  possible.  During  the  course  of  the  cryogenic  programme  it  was  found  that  the  non-linear 
response  of  the  acousto-optic  material  could  cause  significant  power-dependent  depletion  of 
the  acoustic  signal  due  to  the  generation  of  higher  harmonics  [3]. 

In  this  report  the  theory  of  harmonic  generation  due  to  acoustic  non-linearities  is  described. 
The  general  theory,  accurate  to  fourth  order  in  the  non-linearity,  is  derived  in  Section  2. 
Various  simplifications  of  the  general  theory  are  given  in  Section  3,  along  with  example 
calculations.  In  Section  4  the  conclusions  of  the  report  and  its  implications  for  the  construction 
of  wide  bandwidth,  cryogenically  cooled  acousto-optic  signal  processing  systems  are  given.  In 
Appendix  A  the  sufficient  conditions  for  employing  a  linearised  definition  of  the  acoustic 
strain  field  are  derived. 


2  THEORY  OF  ACOUSTIC  NON-LINEARITIES 


i 


t 


The  acoustic  field  propagating  in  the  crj'stal  results  in  a  particle  displacement  denned  by 

U  =  Xo-X«  (1) 


where  xq  is  the  position  of  the  particle  in  its  rest  state  and  Xg  is  the  coordinate  under  the 
influence  of  the  applied  strain  field.  It  has  been  shown  [4]  tha'  for  the  case  of  a  purely 
longitudinal  mode,  the  equation  of  motion  for  the  particle  displacement  is  given  by 


P(r-~ - - 

at2  ax2 


M2  +  Mj> 


ax  '.ax  f 


(2) 


where  po  is  the  density  of  the  crystal  in  the  unstressed  condition  and  x  is  the  Lagrangian 
coordinate  in  the  direction  of  acoustic  propagation.  Equation  (2)  assumes  that  all  harmonics 
1  generated  are  in  the  same  direction  as  the  fundamental  and  are  also  purely  longitudinal 
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modes.  The  Mj  are  non-linear  coefficients,  which  are  combinations  of  elastic  coefficients  up  to 
and  including  order  i  [4] 


Acoustic  measurements  undertaken  by  Signal  Processing  Group  rely  on  exploiting  acousto-optic 
diffraction  of  a  laser  beam  from  the  acoustic  signal  to  obtain  information  about  the  behaviour 
of  the  acoustic  field  [1,  2].  The  intensity  of  the  diffracted  light  is  proportional  to  the 
intensity  of  the  acoustic  strain  field.  As  such,  for  the  purpose  of  comparing  the  final 
theoretical  results  with  experimental  measurements,  it  is  useful  to  convert  Equation  (2)  into 
an  expression  involving  the  one  dimensional  acoustic  strain,  defined  by  [5] 


It  is  shown  in  Appendix  A  that  for  the  case  where  the  non-linearity  is  dominated  by  material 
parameters  satisfying  M2  <  M3  <  M4  <  ...  etc,  it  is  permissible  to  retain  only  the  linear 
contribution  from  Equation  (3).  Substituting  this  into  Equation  (2)  gives 


'S  XI  ^  ..  dlc2dS\ 

—  +  M3 —  S —  +  M4 —  S  — I  + 

3x\  dxl  dx\  dxl 


The  first  term  on  the  right  hand  side  is  just  the  usual  Hooke's  Law  for  one-dimensional 
motion.  The  higher-order  harmonics  arise  as  a  result  of  the  second,  third,  etc  terms  in 
Equation  (4).  Assuming  a  dispersionless  medium,  the  total  strain  field  resulting  from  N 
harmonics  can  be  written 

St  =iy  (Sfxo  +  Spo,e-P^“'-'‘*>)  (5) 

2p=, 

where  Spa  is  the  complex  strain  envelope  at  frequency  po)  and  pk  is  the  corresponding 
wavenumber.  It  is  convenient  to  introduce  the  shorthand  notation 

ST=li(SpePi  +  S‘pe-P’)  (6) 

2p=i 


where  the  term  exp((ot-kx)  is  understood. 

Following  the  procedure  of  Elston  and  Kellman  [5],  equation  (6)  is  substituted  into  Equation  (4) 
and  terms  of  equal  exponential  phase  and  frequency  are  identified.  Taking  the  time 
derivatives  of  Equation  (6)  yields  for  the  left  hand  side  of  Equation  (4) 

=  J-X  ((pw)^  Sp e^  +  c.c.)  (7) 

dt^  2  p., 


The  spectral  component  at  frequency  po)  is  thus 
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= J-(pa))^Sp  ePi 


2  2 


(8) 


a“Sp 

Higher  order  terms  in  — ^  have  been  set  to  zero  on  the  right  hand  side  of  Equation  (7).  In 
ot 

carrying  out  the  spatial  derivatives,  it  is  convenient  to  initially  calculate  only  the  simple 
harmonic  (elastic)  contribution  from  the  first  term  on  the  right  hand  side  of  Equation  (4).  This 
is 


^  =  f^eP'-2kpj^  ePi-(kp)2SpePi  +  c.c. 

3x2  2p.i^3,2  "’’dx 


(9) 


The  spatial  derivatives  are  now  simplified  by  invoking  the  slowly-varying  envelope 
approximation  (SVEA)  i.e. 


k2Sp»k^»^ 

3x2 


(10) 


which  gives  for  the  pco  component. 


32Sp  _ 
3x2  "■ 


(11) 


Substituting  equations  (8)  and  (11)  into  equation  (4)  yields  for  the  normal  elastic  contribution 

-  Po^Ef^Sp  =  -  Sp  +  kpj-^E.|  +  (higher  order  non-linear  terms}  (12^ 

2  \  2  3x  / 


As  previously  noted  the  coefficient  M2  is  the  usual  Hooke's  Law  elastic  coefficient,  and  is 
related  to  co,  po,  k  and  the  acoustic  phase  velocity  Vp  via  the  relations 


Vp=l!i  =  (Ml 
k  *  po 


1/2 


(13) 


Equation  (12)  can  therefore  be  expressed  as 

0S 

— L  =  — —  X  ( higher  order  non-linear  terms  ]  (14) 

pok'pVp' 

Having  exploited  the  simplification  arising  through  the  determination  of  the  harmonic 
contribution,  it  is  appropriate  to  turn  to  the  evaluation  of  the  non-linear  contributions. 
Invoking  the  SVEA  again  yields  expressions  for  the  higher  order  terms  given  by 
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dx  2  m=i 

St  ^  =  - 1  f  X  mjk  (Sm  S„  e*'"*”*  -  S„  +  S„,  Sn*  S„* 

3x  4  ni=l  n=:l 

s?^=-iX  X  rnSnSfC^  ^  -  SniSnSr  ^  +  SmSnSr  ^ 

9x  8  m5=I  n=l  r=l 

-  +  StnSnS;  -  S^SnS;  +  SjuSrSr 

-s;,s;3r‘ 


Differentiating  Equations  (15),  (16)  and  (17)  under  the  constraints  of  the  SVEA,  and 
substituting  the  result  into  equation  (14)  yields 

[Mjj,  X  m  (SmSn  5[p-(m+n)]  -  SmSo  6[p-(n-in)]  +  SmSn  8[p-(m-n)]) 

4poVp  m=l  n=l 

X  Z  S  SmSnSr  5[p-(m+n+r)]  -  SmSnSr  6[p-(n+r-m)] 

2  [11=:  1  n=l  f  =  l  (18) 

+  SmSiSr  5[p-(m+r-n)]  -  SmSi^r  5[p-{r-m-n)]  +  SmSmSr  8[p-(m+n-r)] 

-SnSnSrS[p-(n-m-r)j  +  SinS^Sr  5[p'(m-n-r)])] 


The  delta  functions  express  the  frequency  matching  constraints  on  harmonics  m,  n  and  r,  and 
can  of  course  be  used  to  kill  one  of  the  summations  in  each  term  if  required.  It  is  convenient 

however  to  retain  these  summations  for  the  moment.  These  constraints,  for  the  physically 

»  »  *  •  • 

meaningful  solutions  with  positive  co,  preclude  inclusion  of  terms  such  as  SaiSn  and  SmSnSr. 
The  exponential  frequency  term  exp(pco)  that  must  appear  on  both  sides  of  Equation  (18)  is 
understood. 

In  order  to  clarify  the  underlying  physics  of  Equation  (18),  the  complex  strain  field  Sp  is 
expressed  in  terms  of  the  real  amplitude  and  phase  via 

Sp  H  Sp  exp(jvp)  (19) 

Substituting  Equation  (19)  into  Equation  (18)  and  separating  the  result  into  the  real  and 
imaginary  parts  yields  (dropping  the  caretes  for  simplicity  of  notation) 
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— ^=-Y0)M3X  X  (*‘"(Vn>+Vn-Vp)  S[p-(in+n)]  -  sin(Vn-Vm-Vp)  8[p-(n-m)] 

9x  m=l  n=l 

+  sin(\|/m-Vn-Vp)  8[p-(m-n)]) 

N  N  N 

- — 1  £  ^  m  SmSoSr  {sin(\i/m+VD+VrVp)  5[p-(m+n+r)] 

2  m=l  n=l  r=l 

-  sin(Vn+Vr-Vm-Vp)  8[p-(n+r-m)]  +  sin(\|/m+V|/r-Vn-Vp)  8[p-(m+r-n)] 

-  sin(yrVm-Vn-Vp)  8[p-(r-m-n)]  +  sin(ym+VD-Vr-\i/p)  8[p-(m+n-r)] 

-  sin(yn-Vm-Vr-Vp)  5[p-(n-m-r)]  +  sin(Vm-Vn-Vr-Vp)  8[p-(m-n-r)]) 

where 


YS— L-r.  (21) 

4poVp 

Equation  (20)  describes  the  evolution  of  the  envelope  of  the  strain  field.  The  corresponding 
equation  for  the  imaginary  terms  is 

N  N 

0  =  M3  X  X  f"SmSn  (cos(Vm+Vn-Vp)  S(p-(m+n)]  -  cos(\|/„-\|/m-Vp)  5[p-(n-m)] 

m=l  n=l 

+  cos(Vm-'|/n-Vp)  8[p-(m-n)]) 

+  —  X  XX"'  SmSnSr  (cos(\ym+Vn+'|/r-v|/p)  8[p-(m+n+r)] 

2  ra=l  nsl  r=l  (22) 

-  cos(Vn+V)/rVm-Vp)  5[p-(n+r-m)]  +  co.s(Vm+Vr-Vn-Vp)  6[p-(m+r-n)] 

-  co.s(VrVKm-Vti'Vp)  5[p-(r-m-n)]  +  cos(Vm+Vn-Vr-Vp)  6[p-(m+n-r)] 

-  cos(v|/n-Vm-Vr-Vp)  S[p-(n-m-r)]  +  cos(\|/m-'l/n-Vr-Vp)  8[p-(in-n-r)]) 

Equations  (20)  and  (22)  are  the  final  expressions  for  the  evolution  of  the  magnitude  of  the 
envelope  of  the  strain  field  at  frequency  p©  as  a  function  of  propagation  distance  x,  accurate 
to  fourth  order  in  the  non-linearity.  The  extension  of  the  equations  to  fifth,  sixth  and  higher 
order  non-linearities  is  obvious. 


Equation  (22)  constrains  the  contributing  terms  to  have  phases  that  combine  such  that  their 

sum  is  always  0t  =  ±|k+i-|rt,  k  =  0,  1,  2,  ..  .  This  demonstrates  the  so-called  "phase-locking" 

phenomena  of  the  harmonics  that  has  been  previously  noted[6].  Having  established  all 
possible  terms  that  may  contribute  to  the  evolution  of  the  spectral  component  Sp  of  the  strain 
(as  determined  by  the  delta  functions  and  equation  (20)),  the  phase  constraints  then  select 
those  terms  that,  in  effect,  ensure  conservation  of  energy  of  the  total  system.  The  phase 
locking  implies  that  the  sine  terms  of  Equation  (20)  must  all  be  constrained  to  a  value  of 
either  +1  or  -1.  As  will  be  shown  explicitly  in  the  following  example,  a  choice  of  +1  for  these 
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terms  leads  to  solutions  that  satisfy  conservation  of  acoustic  power  and  are  intuitively  correct. 
Equation  (20)  can  thus  be  written  as 

^  =  -y(0M3  rnSmSn  (S[p-(in+n)]  -  8[p-(n-m)]  +  8[p-(in-n)]) 

m=l  n=l 

-  ^  5^  mSmSnSr(S[p-(in+n+r)]  -  8(p-(n+r-m)]  +  5[p-(m+r-n)]  (23) 

2  m=l  n=l  f=l 

-  S[p-(r-m-n)]  +  2[p-(m+n-r)]  -  5[p-(n-m-r)]  +  8[p-(m-n-r)]) 

In  order  to  demonstrate  conservation  of  acoustic  power  of  the  general  equations  and  the 
validity  of  the  choice  of  sign  for  the  phase  terms.  Equation  (20)  will  be  used  in  the  following 
example. 

Example:  fourth  order  non-linearity,  N=3 
The  contributing  terms  are  found  to  be 

^  =  -ycoM3(SiS:  sin^i  +  S2S3  sin^2)  -  Y5^Mi{sfS3  sin43  +  S2-S3  sin^)  (24) 

3x  2 

sin^i  -  2S1S3  sin^2)  +  2ya)  M4  S1S2S3  sin^4  (25) 

dx 

^  =  3y<bM3  S1S2  sin42  -  I^5Mi/3SiS2  sin^4  -  S?  sin^i)  (26) 

dx  2 

where  the  are  the  combinations  of  phase  terms  appearing  in  Equations  (20)  and  (22). 
Conservation  of  power  requires 


S?  +  S2  +  S3  =  constant 

which  implies  that 

Si^+  S2— +  S3— =  0 
dx  dx  dx 

Substituting  Equations  (24)  to  (26)  into  Equation  (28)  confirms  that  Equations  (24)  to  (26) 
satisfy  conservation  of  acoustic  energy  independent  of  the  choice  of  +1  or  -1  for  the  value  of 
the  sine  terms.  Choosing  the  sine  terms  to  have  a  value  of  +1  is  intuitively  reasonable  since 
(i)  this  leads  to  an  overall  negative  derivative  for  the  fundamental  indicating  that,  at  least 
initially,  the  fundamental  will  suffer  depletion  (as  it  must),  and  (ii)  terms  with  m+n=p  are 
always  associated  with  a  positive  derivative,  as  expected  for  lower  frequency  contributions 
that  "pump "  the  higher  harmonic.  These  considerations  finally  give  a  set  of  equations  for 
the  evolution  of  the  three  harmonics,  accurate  to  fourth  order  in  the  non-linearity,  as 


(27) 
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^  =  -ycoM  3(8182  +  8283)  -  I^^Sf83  +  S283)  (29) 


^  =  yo)M3(8i  -  28183 )  +  2y(oM48i82S3  (30) 

dx 


3YC0M3S182  -  l^^li(3Si82  -  8?) 


(31) 


i 


3  SIMPLIFICATION  OF  THE  GENERAL  THEORY 

In  this  section  several  examples  are  given  for  '.vhich  the  specific  form  of  the  rate  equations 
which  describe  harmonic  generation  due  to  non-linear  acoustic  interactions  can  be  simplified 
under  various  approximations  .  The  examples  given  are 

( i )  the  rate  equations  and  analytic  solution  for  the  third  order,  two  harmonic  case 

(ii)  a  simplified  form  of  the  rate  equations  for  the  third  order,  N  harmonic  case 

( i  i  i )  parameterisation  of  the  full  solution  for  low  to  moderate  acoustic  powers 

Acoustic  attenuation  is  included  in  a  phenomenological  manner  in  each  of  the  cases  as  a 
simple  additional  term  in  the  rate  equations.  For  the  purpose  of  the  calculations  the  L[001] 
propagation  mode  of  Te02  at  a  frequency  of  1  GHz  is  used  as  the  standard  example,  and  an  ap- 
frequency  dependence  is  assumed.  A  transducer  with  radius  0.125  mm  is  used  throughout  and 
the  standard  input  power  for  all  calculations  is  15  mW  unless  otherwise  specified. 

The  L[(X)1]  mode  of  Te02  exhibits  an  anomalously  large  acoustic  non-linearity  due  to  a  weak 
oxygen-oxygen  bond  in  the  lattice  [7],  and  is  therefore  ideally  suited  to  the  study  of  acoustic 
non-linear  phenomena.  Furthermore,  despite  the  fact  that  Te02  is  piezoelectric,  the  L[001] 
propagation  mode  is  not  effected  (8]  and  hence  there  is  no  piezoelectric  stiffening  of  the 
elastic  constants  of  the  material. 

Numerical  solution  of  the  coupled  differential  equations  that  occur  in  examples  (i)  and  (ii) 
was  undertaken  using  a  fourth  order  Runge-Kutta  method. 

( i )  Third  order,  two  harmonic  case: 

In  this  case  the  rate  equations  simplify  to 

^=-aiSi  -M3YO)SiS2  (32) 

^x 
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^  =  -0282  +  MjymS?  (33) 

bx 

where  a]  2  denotes  the  acoustic  absorption  coefficient  in  Nepers/m  for  frequency  o)]  2- 
Equations  (32)  and  (33)  can  be  solved  analytically  for  the  case  ai  2  =  0,  resulting  in 


Si(x)  =  S!{0)  sech(Y(oSi{0)Mix) 

(34) 

S2(x)  =  Si{0)  tanh(Ya)Si(0)M3x) 

(35) 

The  strain  is  related  to  acoustic  power  Pg  by 

=  (36) 

2 

where  A  is  the  cross  sectional  area  of  the  acoustic  beam.  Thus,  Equations  (34)  and  (35)  can  be 
written  in  terms  of  the  ii  itial  acoustic  power  in  the  fundamental,  Pq,  as 

Pi(x)  =  Posech^ryPo/A  cox)  (37) 

P::(x)  =  Potanh2(rVPo/AcDx)  (38) 

where  the  harmonic  generation  constant  [5],  T,  is  defined  by 

r  = - (39) 

(2Po)‘’V^^ 

Figure  1  demonstrates  second  harmonic  generation  due  to  the  third  order  non-linearity  of  the 
L[CX)1]  mode  of  Te02,  for  a  frequency  of  1  GHz.  From  Equations  (37),  (38)  and  (39)  it  is 

apparent  that  non-linear  effects  will  be  important  for  Bragg  cells  that  operate  at  high 
frequencies  and  have  slow  phase  velocities  (and  hence  large  time-bandwidth 
products)  -  unfortunately  the  very  properties  that  are  highly  desirable  for  EW/SIGINT 
applications.  Figure  2  shows  the  numerical  solution  of  Equations  (32)  and  (33)  for  die  same 
parameters  as  in  Figure  1  but  including  an  attenuation  of  0.5  dB/mm,  assuming  an  co^  frequency 
dependence  for  the  attenuation  coefficient.  The  effect  of  attenuation  is  clearly  to  decrease  the 
coupling  of  the  fundamental  h  he  harmonic,  and  extend  the  distance  over  which  significant 
power  can  be  found  in  the  fundamental.  This  suggests  that  in  systems  dominated  by  acoustic 
nonlinearities  a  small  amount  of  frequency  dependent  attenuation  is  in  fact  beneficial  to  the 
propagation  characteristics  of  the  fundamental. 

( i  i )  Third  order  non-linear..y,  .N  harmonic  case; 

Probably  the  easiest  method  to  obtain  the  general  form  of  the  rate  equations  for  an  arbitrary 
number  of  frequencies  is  that  of  writing  down  the  equations  for  the  first  M  harmonics  and 
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generalising  from  the  observed  symmetry  of  the  resulting  coupled  differential  equations.  For 
example,  taking  the  fundamental  and  first  five  harmonics  yields 


=  -otiSi  -  y{i3M3(SiS2  +  S2S3  +  S3S4  +  S4S5  +  S5S6)  (40) 

dx 

^  =  -a2S2  +  2Ya)M3|^  -  S1S3  -  8384  -  8385  -  8486)  (4 1 ) 

dx  \  2  I 

^  =  -a3S3  +  3Ya)M3(  S1S2  -  8184  -  S2S5  -  S3S6)  (42) 

dx 

^  =  -0484  +  4Ya)M3(  S183  +  ^  -  S185  -  82S6j  (43) 

^  =  -asSs  +  5yo)M3(  S1S4  +  S2S3  -  8,85)  (44) 

dx 

^  =  -0686  +  6Ya)M3|  S1S5  +  S2S4  +  (45) 

dx  '  2  f 


Ignoring  the  loss  terms,  the  symmetry  of  these  equations  clearly  demonstrates  the  following 
rules  for  writing  down  the  rate  equation  describing  the  evolution  of  the  harmonic  at  frequency 
pcd.  These  are,  for  each  term  SmSn(noting  that  m  <  n  always) 

( 1 )  if  n  <  p  then  the  term  is  +5^3^ 

(2)  if  n  >  p  then  the  term  is  -S^Sn 

q2 

(3)  if  m  =  n  then  the  term  is  +  ^ 

2 

These  observations  lead  to  the  general  third  order  expression  for  the  harmonic  at  frequency 
pco  as 

(  L  N-p  \ 

^  =  -OpSp  +  pywmJ  (1-6, p)  X  SmSNLBL  .  (i-S^p)  ^  (46) 

dx  \  (ji_]  1  +  [n=l  / 

where  the  5ij  are  Kronecker  deltas,  t  =  £  (only  integer  for  p  even)  and  L  =  TRUNC(t),  where 
TRUNC  means  "truncate  to  integer". 

Figure  3  shows  the  behaviour  of  the  fundamental  calculated  from  the  numerical  solution  for 
the  L(0011  mode  for  the  cases  of  N  =  2, 10,20  and  40  harmonics  and  no  attenuation.  It  is 
apparent  that  even  in  the  relatively  short  space  of  a  typical  Bragg  cell  aperture,  the  effect  of 
the  very  high  harmonics  on  the  behaviour  of  the  fundamental  is  substantial,  indicating  that 
under  low  loss  conditions  accurate  modelling  of  the  fundamental  becomes  very  involved. 


UNCLASSIFIED 


9 


ERL-0643-RR 


UNCLASSIFIED 


Another  interesting  feature  of  Figure  3  is  that  the  energy  does  not  completely  drain  from  the 
fundamental  into  higher  harmonics,  as  is  suggested  by  the  results  for  N=2.  In  fact,  as  shown  in 
Figure  4,  after  a  sufficiently  large  propagation  distance  a  dynamic  equilibrium  of  the 
distribution  of  acoustic  power  between  the  N  modes  is  achieved.  Figure  4  shows  the  full 
numerical  solution  of  Equations  (40)  to  (45)  for  the  case  of  zero  attenuation. 

The  effect  of  an  acoustic  attenuation  of  0.(X)2  dB/mm  is  demonstrated  in  Figure  5.  As  expected 
the  assumed  afi  dependence  of  the  attenuation  very  significantly  reduces  the  influence  of 
higher  harmonics,  and  it  is  interesting  to  note  that  the  oscillations  in  the  fundamental  are  not 
only  damped  but  the  periodicity  is  also  slightly  altered  in  the  presence  of  finite  attenuation. 

It  is  clear  from  these  results  that  any  acousto-optic  signal  processing  device  relying  on  a 
cryogenically  cooled  Te02  Bragg  cell  operating  in  the  GHz  regime  would  have  its 
performance  severely  effected  by  acoustic  non-linearities.  This  is  expected  to  also  be  true  of 
materials  such  as  GaP  and  LiNbOs,  which  could  begin  to  exhibit  similar  behaviour  at  around 
the  10  GHz  region  (the  non-linear  coefficients  of  these  materials  are  about  an  order  of 
magnitude  smaller  than  those  of  the  L[001]  mode  of  Te02  [5],  and  as  such  will  only  play  a  role 
at  higher  frequencies.  At  such  high  frequencies  however,  dispersion  effects  could 
significantly  modify  the  non-linear  response). 

This  degradation  in  the  high  frequency  performance  of  acousto-optic  devices  is  intrinsic  to  the 
physical  nature  of  the  material  and  as  such,  unlike  the  deleterious  effects  of  thermal 
attenuation,  represents  a  fundamental  limitation  to  the  performance  of  these  devices  in  the 
GHz  regime. 

One  possibility  for  operating  acousto-optic  devices  at  multi-gigahertz  frequencies  as 
suggested  by  Figures  (4)  and  (5)  may  be  to  exploit  the  region  in  which  there  is  a  dynamic 
equilibrium  between  the  harmonics.  This  could  be  achieved  by  constructing  longer  cells  (an 
option  made  feasible  by  the  negligible  attenuation  at  cryogenic  temperatures),  or  by  operating 
at  high  frequencies  or  high  powers  (the  latter  options  essentially  "compress"  the  highly 
oscillatory  region).  The  effects  of  acoustic  non-linearities  may  then  be  no  more  troublesome 
than  the  more  familiar  near-field  integrated  optical  effect  [9],  which  causes  anomalous 
diffraction  of  light  in  regions  very  near  the  transducer.  Thus,  for  example,  by  exploiting  the 
region  of  dynamic  equilibrium  it  may  be  possible  to  construct  a  very  high  frequency  acousto¬ 
optic  deflector  suitable  for  mode-locking  diode  lasers  or  for  noise  suppression  applications.  In 
most  signal  processing  applications  however,  depletion  of  the  fundamental  due  to  the  non¬ 
linear  response  of  the  acousto-optic  material  will  fundamentally  limit  the  performance  of 
ciy'ogenically  cooled,  high  frequency  devices. 

( i  i  i )  Parameterised  model  for  the  non-linear  behaviour: 

In  section  3(i)  it  was  shown  that  the  most  fundamental  non-linear  interaction  involves  a  sech^ 
depletion  of  the  fundamental  to  the  second  harmonic  (Equation  (37)).  Thus,  for  low  powers,  it 
is  reasonable  to  assume  that  the  power  depletion  for  the  general  case  should  initially  be 
dominated  by  this  sech^  behaviour.  From  Figure  (4)  it  is  apparent  that  the  major  effect  of 
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accounting  for  higher  harmonics  in  the  model  is  to  allow  acoustic  energy  to  return  to  the 
fundamental,  resulting  in  a  slowing  of  the  overall  rate  of  depletion  and  a  finite  amount  of 
energy  always  being  present  in  the  fundamental.  From  Figure  (5),  it  can  be  seen  that  the  effect 
of  acoustic  attenuation  is  to  reduce  this  feedback  and  finally  to  dominate  the  decay 
characteristics  for  large  attenuations  or  large  propagation  distances.  Such  observations 
suggest  a  very  simple  model  for  the  decay  of  power  in  the  fundamental,  ie 

Pi(x)  =  Posech^rVPo/A  0)x)  +  Pb  (47) 

where  the  addition  of  the  baseline  Pg  is  the  simplest  method  of  retaining  energy  in  the 
fundamental.  Thus,  the  parametric  model  basically  describes  a  sech^  decay  of  the 
fundamental  with  the  baseline  (for  which  there  is  no  justification  in  the  exact  theory) 
accounting  for  the  influence  of  coupling  to  higher  harmonics. 

The  result  of  fitting  Equation  (47)  to  the  full  theoretical  curves  is  shown  in 
Figures  (6)  and  (7).  An  attenuation  of  0.004  dB/mm  (which  corresponds  roughly  to  a 
temperature  of  4  K)  has  been  used  in  the  calculation  of  the  exact  solutions.  Figure  (8)  shows  a 
plot  of  the  acoustic  power  based  on  the  parametric  model  as  a  function  of  acoustic  power  used 
in  the  exact  solutions.  It  is  clear  from  these  results  that  the  parametric  model  describes  the 
actual  behaviour  quite  well  for  low  to  moderate  acoustic  powers  using  the  indicated  scaling 
parameters.  This  analytic  model  can  therefore  be  used  to  describe  the  depletion  of  the 
fundamental  for  low  to  moderate  acoustic  powers  and  small  attenuation,  and  provides  a 
simple  model  for  the  behaviour  of  the  fundamental  in  the  presence  of  acoustic  non-linearities. 
Such  a  model  can  be  applied  advantageously  in  characterising  experimental  data,  where 
only  the  parameters  T  and  Pq/Pb  need  be  recorded. 


4  CONCLUSIONS 

In  this  report  an  extension  of  the  third  order  theory  of  acoustic  non-linearities  based  on  the 
approach  by  Elston  and  Kellman  [5]  has  been  derived.  The  necessity  of  going  to  fourth  order 
in  the  non-linear  response  was  driven  by  the  requirement  for  a  theoretical  model  capable  of 
describing  the  experimentally  observed  characteristics  of  the  power  dependent  depletion  of  a 
signal  propagating  in  a  cryogenically  cooled  acousto-optic  material  exhibiting  extremely  low 
attenuation.  These  results  have  been  reported  elsewhere  [3]. 

In  Section  2  the  general  expression  for  the  strain  field  accurate  to  fourth  order  in  the  non¬ 
linearity  was  derived. 

In  Section  3  several  specific  examples  were  given  which  highlighted  simplified  descriptions 
for  tlie  non-linear  response,  and  a  parametric  model  accurate  for  low  to  moderate  acoustic 
powers  was  presented. 

The  results  of  the  modelling  indicate  that  power  dependent  depletion  of  the  input  signal  due 
to  harmonic  generation  will  be  a  fundamental  limitation  to  the  utility  of  acousto-optic  signal 
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processing  architectures  operating  in  the  gigahertz  regime.  Acousto-optic  materials  with 
slow  phase  velocities  capable  of  operating  at  high  frequencies  (i.e.  those  materials  that 
have  a  large  time-bandwidth  product,  which  corresponds  to  high  resolution)  are  the  most 
severely  affected. 

Under  certain  circumstances,  the  region  of  dynamic  equilibrium  which  is  established  between 
harmonics  may  find  some  use,  however  it  is  unlikely  that  these  will  include  the  very  high 
resolution  requirements  of  many  EW/SIGINT  signal  processing  applications. 
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Figure  1  Analytic  solution  of  the  two  harmonic,  third  order  non-linear  depletion  of  the 

fundamental  for  the  L[00l]  mode  of  Te02. 
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Figure  6  Parameterised  fit  to  full  numerical  solution  for  low  acoustic  input  power. 
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Figure  7 


Parameterised  fit  to  full  numerical  solution  for  higher  acoustic  input  power. 
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Figure  8  Relaiionship  between  the  analysed  input  power  Pq  for  the  parameterised  model  as 
a  functic  a  of  P©  from  the  full  numerical  solution. 
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APPENDIX  A 

JUSTinCATION  OF  THE  LINEAR  STRAIN  APPROXIMATION 


In  this  Appendix  the  condition  under  which  it  is  valid  to  retain  only  the  linear  contribution 
to  the  strain  field  is  derived.  The  starting  point  is  the  non-linear  equation  of  motion  for  the 
particle  displacement 


Po- 


^^u  _  9^u 


M2  +  Ms  —  +  M4  ( — f  +  Ms  ( — f 

3x  \dxl  \dxl 


(Al) 


As  indicated  in  Section  2,  the  strain  field  for  one-dimensional  propagation  is 

dx  2 \3x  I 

Defining 

=  ~  thenC  =  S-^ 

ax  2 


i.e. 


Then  using  =  2(S  -  Q  on  the  L.H.S.  gives 

^^s-is^+is^  +  oiC*) 

2  2 


(A2) 


(Al) 


The  term  in  square  brackets  in  Equation  (Al)  is  then  of  a  form,  accurate  to  third  order  in  the 
strain,  given  by 
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Now, 


M2  +  Nt;  +  =  M2  +  M?  js  -  W  j  +  Mt  (s^  -  s^) 

=  M2  +  M3S  +  jM4-^js^  +  j^-M4 


^=ii=^-s^  +  2s  — 

3x  3x  dx  2  3x 

M2  +  M3^4-Mi— 1  +  ••• 

dx  \3x| 

=  j^.S^  +  3s2^Vi^  +  M3S  +  (M4-M3js^  +  \ 

\^x  dx  2  dx^  \  2  I  I 

=  M2^  +(M3  -  M2)S—  +  (3.M2  -  2-M3  +  Ml  )s^—  +  o(s^). 

dx  dx  h  2  r  ax 


.  ^ 

"dx^ 


(A4) 


(A5) 


Substituting  Equation  (A5)  into  the  R.H.S.  of  Equation  (Al)  and  differentiating  with  respect 
to  X  yields 


au 

Po — r 

at2 


^  =  M2^  +  (Ms  -  M2)— (s— 1  +  (2.M2  -  2.M3  +  M4) — (s2— j  ( A6) 

ax  3^2  '  dx\  3x1  \2  2  ^3x1  dxf 


The  term  on  the  left  hand  side  of  Equation  (A6)  is,  from  Equation  (A3),  of  the  form 


(A7) 


Only  the  linear  term  in  S  need  be  retained  since  typically  S-IO'^,  and  as  such  the  higher 
order  terms  are  very  small.  Furthermore,  if  the  material  constants  satisfy 


M2  «  Ms  «  M4  «  .  .  . 


(A8) 


(as  is  usually  the  case  if  acoustic  non-linearities  are  to  be  seen  at  all,  due  to  the  very  small 
magnitude  of  the  strain),  then  Equation  (A6)  can  be  expressed  as 


Po-^  =  M; 

3t^ 


3^S 


dx 


+  Ms— 
2  3x 


\  dxj  3x  \  dx I 


{A9) 


Equation  (A9)  is  identical  to  Equation  (4)  of  Section  2,  which  is  obtained  directly  from  the 
approximation  of  using  only  the  linear  contribution  to  the  strain  field.  The  difference 
between  Equation  (4)  and  Equation  (A9)  is  that  the  strain  in  Equation  (A9)  is  the  full,  non¬ 
linear  strain  defined  by  Equation  (A2).  It  is  clear  that  the  material  constraints  defined  by 
Equation  (A8)  are  a  sufficient  condition  to  allow  use  of  the  linearised  expression  for  the  strain 
field. 


24 


UNCLASSIFIED 


UNCLASSIFIED 


ERL-0643-RR 


DISTRIBUTION 

Copy  No. 

Defence  Science  and  Technology  Organisation 

Chief  Defence  Scientist  ) 

Central  Office  Executive  )  1 

Counsellor,  Defence  Science,  London  Cont  Sht 

Counsellor,  Defence  Science,  Washington  Cont  Sht 

Scientific  Adviser,  Defence  Central  1 

Scientific  Adviser,  Defence  Intelligence  Organisation  1 

Navy  Scientific  Adviser  1 

Air  Force  Scientific  Adviser  1 

Scientific  Adviser,  Army  1 

Electronics  Research  Laboratory 

Director  1 

Chief,  Information  Technology  Division  Cont  Sht 

Chief,  Electronic  Warfare  Division  Cont  Sht 

Research  Leader,  Signal  and  Information  Processing,  EWD  1 

Dr.  M.L.  Brown,  EWD  1 

Chief,  Communications  Division  1 

Research  Leader,  Communications  Countermeasures  1 

Dr  A.C.  Lindsay  1 

Mr  S.C.  Troedson 

Media  Services  1 

Libraries  and  Information  Services 

Australian  Government  Publishing  Service  1 

Defence  Central  Library,  Technical  Reports  Centre  1 

Manager,  Document  Exchange  Centre,  (for  retention)  1 

National  Technical  Information  Service,  United  States  2 

Defence  Research  Information  Centre,  United  Kingdom  2 

Director  Scientific  Information  Services,  Canada  1 

Ministry  of  Defence,  New  Zealand  1 

National  Library  of  Australia  1 

Defence  Science  and  Technology  Organisation  Salisbury,  Research  Library  2 

Library  Defence  Signals  Directorate,  Melbourne  1 

British  Library  Document  Supply  Centre  1 

Spares 

Defence  Science  and  Technology  Organisation  Salisbury,  Research  Library  6 


UNCLASSIFIED 


25 


ERL-0643-RR 


UNCLASSIFIED 


26 


UNCLASSIFIED 


Department  of  Defence 

DOCUMENT  CONTROL  DATA  SHEET 

Page  Classification 

UNCLASSIFIED 

Privacy  Marking/Caveat 
(of  Document) 

N/A 

la.  AR  Number 
AR-006-986 

lb.  Establishment  Number 

ERL-0643-RR 

2.  Document  Date 

AUGUST  1992 

3.  Task  Number 

DEF  89/210 

4.  Title 

AN  ANALYSIS  OF  THE  EFFECT  OF 

5.  Security  Classification 

6.  No.  of  Pages 

26 

[Zl  IZ]  in 

7.  No.  of  Refs. 

9 

WIDE  BANDWIDTH  ACOUSTO-OPTIC 
BRAGG  CELLS 


S  (Secret)  C  (Conti )  R  (Rest)  U  (Unclass) 

*  For  UNCLASSIFIED  docs  with  a  secondary  distribution 
LIMITATION,  use  (L)  in  document  box. 


8.  Author(s) 


Anthony  C.  Lindsay,  Ian  G.  Fuss  and 
Shaun  C.  Troedson 


9.  Downgrading/Delimiting  Instructions 

N/A 


10a.  Corporate  Author  and  Address 

Electronics  Research  Laboratory 
PO  Box  1 500 
SALISBURY  SA  5108 


10b.  Task  Sponsor 


1 1 .  Officer/Position  responsible  for 

Security . 

Downgrading . 

Approval  for  Release . D.ERU... 


DEF 


1 2.  Secondary  Distribution  of  this  Document 

APPROVED  FOR  PUBLIC  RELEASE 

Any  enquiries  outside  stated  limitations  should  be  referred  through  DSTIC,  Defence  Information  Services, 


^  - ‘ -  ■ - ■ - — - 

i 

13a.  Deliberate  Announcement 

No  limitation 

j  1 3b.  Casual  Announcement  (tor  citation  in  other  documents) 

i 

1  ^  1 

No  Limitatbn 

Ref.  by  Author ,  Doc  No.  and  date  only. 

14.  DEFTEST  Descriptors 

Signal  processing.  Acousto-optics,  Acoustic  non-linearities* 

15.  DISCAT  Subject  Codes 

1711 

16.  Abstract 


In  this  report  an  extension  of  the  theory  of  acoustic  non-linearities  accurate  to  fourth 
order  is  derived.  A  new  fundamental  limitation  on  the  high  frequency  operation  of 
acousto-optic  signal  processing  architectures  is  demonstrated. 


Page  Classification 

UNCLASSIFIED 


17.  Imprint 


Electronics  Research  Laboratory 
PO  Box  1500 
SALISBURY  SA  5108 


1 8.  Document  Series  and  Number 


19.  Cost  Code 


20.  Type  of  Report  and  Period  Covered 


ERL-0643-RR 


723238 


RESEARCH  REPORT 


22.  Establishment  File  Reference(s) 


23.  Additional  informatbn  (if  required) 


Doc  S9Cf 


